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ABSTRACT 


» 

We  examine  the  validity  of  the  application  of  weak  scatter  theory  to 
t  ionospheric  scattering  data.  We  conclude  that  the  weak  scatter  theory  is 

not  universally  applicable.  We  also  compare  Born  and  Rytov  weak  scatter 
calculations  and  find  the  latter  superior.  The  Rytov  calculations  also 
I  appear  to  explain  the  enhancement  in  the  scintillation  index  S^/c^  in 

the  equatorial  vis  a  vis  the  polar  stations. 
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Random  fluctuations  or  scintillations  of  the  phase  and  amplitude  of 

electromagnetic  signals  propagating  through  the  Ionosphere  have  been  observed 

experimentally.  These  signal  fluctuations  are  attributed  primarily  to  the 

fluctuations  In  the  plasma  density  In  the  F- layer  of  the  Ionosphere  and  depend 

on  such  factors  as  the  wavelength  of  the  signal,  the  source-receiver  geometry, 

the  geo-physical  conditions,  etc.  In  the  past  two  decades  both  in  situ  and 

remote  observations  have  been  made  to  elucidate  the  morphological  properties 

of  irregularities  causing  the  scintillation.  The  earliest  observations  of 

the  scintillations  Involved  radio  emissions  from  stellar  sources^.  These 

2 

scintillations  were  later  Identified  as  originating  In  the  ionosphere  . 

Booker3  has  made  a  detailed  review  of  the  radio  star  observations.  With  the 

advent  of  satellites,  controlled  experiments  could  be  performed.  Yeh  and 

Swenson4  have  reviewed  the  early  satellite  measurements.  In  the  mid-  and  late 

seventies  wide  band  satellites  In  polar  orbits  were  launched  and  a  vast 

5-13 

amount  of  detailed  Information  gathered  .  Particularly  relevant  to  this 
report  are  the  data  obtained  from  the  U.S.  Navy  Navigational  Systems  satel¬ 
lites®’^  and  the  Air  Force  P76-5  wideband  satelllte^”^. 

Theoretical  analysis  of  the  scintillation  problem  has  been  carried  out 

14 

by  several  authors.  Booker  et  al.  proposed  that  the  diffraction  of  the 
radio  waves  by  the  fluctuations  In  electron  density  In  the  Ionosphere  could 
cause  the  observed  scintillation.  They  modeled  the  effect  of  the  electron 
density  fluctuations  by  a  thin  phase  changing  diffraction  screen.  Further 
work  on  the  phase  screen  model  was  carried  out  by  Mercler^5  and  Briggs  and 
Parkin^6.  The  scintillation  problem  has  also  been  analyzed  using  the  weak 
scintillation  models  developed  by  researchers  In  other  fields  (see,  for 
example,  the  monograph  of  Tatarski  ^).  Under  weak  scintillation,  either  the 
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Born  approximation  or  the  Rytov  approximation  solution  to  the  scalar  wave 

equation  Is  applicable.  Crane  6,7  and  Rlno12  and  Fremouw13  have  used  such 

weak  scatter  formulas  to  analyze  the  satellite  data.  Strong  scatter  theories 

18  19  20 

have  been  developed,  among  others,  by  Prokhorov  et  al.  ,  Fante  ,  Ishimaru 

21  22  23 

and  Rumsey  .  Recently,  Dashen  *  has  developed  a  path  integral  formulation 
of  the  problem. 

These  above  mentioned  approaches  can  be  used  to  yield  such  measures  of 
scintillation  as,  the  variance  of  phase  ^$2^,+  the  variance  of  amplitude^2/* 
variance  of  intensity  S^,  scintillation  index  etc.  These  are  basically 

various  types  of  moments  of  the  distribution  function  of  the  electric  field. 

To  date,  there  has  been  no  satisfactory  prediction  of  the  distribution  of  the 
field  for  a  general  ionospheric  scattering  situation.  For  a  weak  scintilla¬ 
tion,  an  argument  based  on  the  central  limit  theorem  has  been  advanced  which 
yields  Gaussian  distributions  for  the  in-phase  and  quadrature  components  of 
the  electric  field  in  the  Born  approximation  and  for  the  phase  and  log- 
amplitude  in  the  Rytov  approximation17.  Crane6*7  has  argued  that  the  use  of 
the  central  limit  theorem  cannot  be  justified  for  what  amounts  to  single 
scattering  approximations.  In  any  event,  the  log-normal  distribution  does  not 
appear  to  fit  the  ionospheric  data  well.  For  practical  data  analysis  the 
Nakagaml-m  distribution  appears  to  be  useful1®’26.  The  Nakagami-m  model  is 
restricted  to  the  description  of  intensity  statistics  and  provides  an  Intensity 
distribution  function  which  reduces  In  the  appropriate  limits  to  the  Rice 
Distribution  and  the  Rayellgh  Distribution,  respectively.  Fremouw  et  al.10 
have  recently  proposed  a  two-component  model  in  which  the  electric  field  Is 

t  The  brackets  denote  an  average  over  the  plasma  fluctuations. 
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modeled  as  the  product  of  a  dlffractively  scattered  component  E.=X  +  iYe  and 
a  refractlvely  focused  component  Ef=exp[Xf  +  1*f].  Xs>  Ys>  Xf  and  are 
assumed  to  be  Gaussian  and  Es  Is  uncorrelated  with  E^.  This  model  requires  six 
Independent  parameters,  which  must  be  found  from  the  data,  for  a  complete 
description. 

The  purpose  of  this  final  report  is  to  examine  the  compatibility  of  the 

ionospheric  data  with  the  various  theoretical  approximations.  The  ionospheric 

scattering  Is  quite  complicated  because  of  the  anisotropy  of  the  irregularities, 

the  power-law  spectrum  of  scales  in  the  plasma  fluctionas  and  the  source  motion. 

These  factors, coupled  with  the  data  processing  constraints, prevent  a  single 

viable  theoretical  approach  for  the  entire  satellite  pass.  Unlike  the  conclu- 
12 

sion  of  Rlno  ,  In  most  cases,  the  scattering  cannot  be  determined  to  be  weak 
for  an  entire  pass.  Most  data  show  Intense  scintillations  for  at  least  a 
portion  of  the  pass,  particularly  in  the  auroral  region.  Hence,  a  strong 
scatter  theory  has  to  be  used  In  these  regions.  Also,  in  agreement  with  Crane's 
analysis®’^,  the  Rytov  approximation  is  preferable  to  the  Born  approximation, 
particularly  In  the  equatorial  region. 

In  Section  2,  we  briefly  review  the  literature  on  the  ionospheric 
scattering.  Section  3  contains  an  analysis  of  the  Rytov  approximation.  In 
Section  4,  the  Rytov  and  Bom  results  for  a  number  of  satellite  pass  configura¬ 
tions  are  compared  and  are  shown  to  be  significantly  different.  Complete 
passes  both  for  equatorial  and  polar  stations  are  calculated  in  the  Rytov 
approximation  and  the  validity  of  this  approximation  is  discussed.  Section  5 
reiterates  some  of  the  Important  conclusions. 
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2. 


A  REVIEW  OF  IONOSPHERIC  SCINTILLATION  THEORIES  AS  APPLIED  TO  THE 
SATELLITE  DATA 


(i)  Helmholtz  and  Parabolic  Equations 

Ionospheric  scintillation  Is  caused  by  electron  density  irregularities 
which  act  at  fluctuations  in  the  refractive  index  for  the  electromagnetic 
wave14.  A  complete  theory  for  wave  propagation  through  a  medium  with  random 
variations  in  the  refractive  Index  is  not  available  but  approximate  solutions 
can  be  sought  in  the  various  parameter  domains.  Development  of  these  approxi¬ 
mate  solutions  has  been  initiated  in  such  related  fields  as  optical  propaga- 

1 7  23 

tion  through  the  atmosphere  ,  acoustical  propagation  in  the  ocean  ,  radiowave 

transmission  through  the  Ionosphere,  etc.  The  work  prior  to  1970  Is  reviewed 

in  the  article  of  Barbanenkov  et  al.2®.  The  work  carried  out  in  the  70' s  is 

19  18  7 

reviewed  by  such  authors  as  Fante  ,  Prokhorov  et  al.  ,  and  Crane  .  The  path 

integral  formulation  of  the  signal  propagation  in  random  media  is  reviewed  in 

a  book  by  Flattd  et  al.23.  Applications  of  these  theories  specifically  to  the 

6  7  *  9-13 

inospherfc  scattering  have  been  carried  out  by  Crane  ’  ,  Fremouw  et  al.  , 

Rlno  et  al.^”13  and  Seshadrl  and  West27. 

The  usual  starting  point  for  the  analysis  of  the  effects  of  a  random 
refractive  Index  on  the  propagation  of  scalar  waves  is  the  scalar  wave  equation. 
The  use  of  the  scalar  wave  equation  to  describe  the  propagation  of  electro¬ 
magnetic  waves  in  the  presence  of  a  wide  range  of  refractive  index  fluctuations 
cannot  be  justified  a  priori.  However,  the  Investigation  of  the  polarization 
effects  in  the  vector  Helmholtz  equation  indicates  that  for  Irregularities 
In  the  far  field  of  both  the  source  and  the  receiver  and  possessing  scales  L»X 

the  wavelength  of  propagation,  the  depolarization  effect  is  second  order  and  may 
28  29  30 

be  neglected  ’.  Booker  and  Ferguson  have  Inferred  that  the  Inner  scale 
cut-off  is  approximately  the  ion  gyro  radius,  ^5m.  From  these  considerations. 
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the  scalar  wave  equation  should  prove  adequate  for  the  largest  wavelengths  of 
Interest  (of  the  order  of  3m),  and  the  polarization  effects  should  prove  negllgibl 
The  scalar  wave  equation  for  the  wave  amplitude  f(r,z,t)  is  given  by 

7%(r,z,t)-V^<r*z»t)  *  0  (2.1) 

c2  at2  “ 

where  £  *  (x,y)  are  the  transverse  coordinates  and  z  is  in  the  direction  of 
propagation.  Assuming  a  solution  of  the  form 

f(_r *z»t)  =  Re{E(r,z,t)  ei(kz_u>t)|  (2.2) 

and  substituting  (2.2)  into  (2.1)  one  obtains  the  Helmholtz  equation 

72E(r,z,t)  +  k2n2(r,z,t)E(r,z,t)  «  0  .  (2.3) 

The  quantity  n(r,z,t)  Is  the  refractive  index  of  the  medium,  which  in  the  case 
of  the  ionosphere.  Is  a  random  function  of  position  due  to  the  presence  of 
electron  density  irregularities. 

The  Helmholtz  equation  (2.3)  Is  an  elliptic  partial  differential  equa¬ 
tion,  so  that  In  order  to  find  the  electric  field  at  a  given  spatial  point  one 
must  solve  the  equation  for  the  field  at  all  points  in  space.  This  property 
limits  the  analytical  and  numerical  tractabllity  of  (2.3).  In  the  weak 
scintillation  case,  perturbative  solutions  such  as  the  Born  and  Rytov  approxi¬ 
mations  can  be  obtained  from  the  Helmholtz  equation.  In  the  strong  scintilla¬ 
tion  case  the  Helmholtz  equation  is  first  approximated  by  the  parabolic  equa¬ 
tion;  (we  will  presently  give  a  mathematical  definition  of  weak  and  strong 
scintillation) 


2M  E(r,z,t|r0rfl)+j7»  *  2k![n!(£.2.t)-l]J  E(_r,z,t| j^.0)  »  0,  (2-4) 

,,i’ere  n  ’  |£ + fjr- 


This  equation  is  the  starting  point  of  the  moment  analysis  of  such  authors  as 

Barbenenkov  et  al.2®,  Tatarski1^,  Taylor31 , de  Wolf32  and  Rumsey21.  It  is  also 

22  23 

the  starting  equation  for  Dashen's  path  integral  formulation  '.  In  the 
parabolic  approximation,  the  normal  to  the  wave  front  of  the  electric  field 
is  assumed  to  remain  close  to  the  z-axis  and  the  field  is  assumed  to  obey  the 
boundary  condition 

lim  E(r,z|^,0)  eik(^}  /2z  .  (2'5) 

2 

The  square  of  the  refraction  index  n  (r,z,t)  is  a  product  of  the  dielec¬ 
tric  coefficient  e(r,z,t)  and  the  magnetic  permeability  of  the  medium.  When 
the  m.  urn  is  turbulent,  as  in  the  case  of  the  ionosphere,  the  refractive  index 
Is  a  fluctuating  quantity.  The  magnetic  permeability  is  usually  assumed  to  be 
constant  and  set  equal  to  its  vacuo  value  of  unity.  The  dielectric  coefficient 
e(x,t),  a  random  function  of  x  3  { r,z )  and  t.  Is  defined  in  terms  of  Its  statis¬ 
tical  properties.  The  average  value  of  e(x,t) ,  where  the  average  Is  performed 
over  an  ensemble  of  realizations  of  the  plasma  fluctuations,  can,  without  loss 
of  generality,  be  set  equal  to  unity.  Thus,  one  can  partition  e(x,t)  as  a  sum 
of  Its  average  value  and  the  fluctuations  u(x,t)  about  this  average,  i.e., 

e(x,t)  ■  1  +  u(x,t)  (2.6) 

where,  by  definition.  Indicating  the  average  by  a  bracket,  we  have 

<u(x,t)>  *  0.  (2.7) 
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The  fluctuations  p(x,t)  are  further  specified  by  their  two-point  correlation 
function  p  as 


p(x,x',t,t')  =  <y(x,t)  y(x',t')>  . 


(2.8) 


for  a  stationary  homogeneous  medium,  P  is  a  function  only  of  the  spatial  dif¬ 
ference  x-X"  and  time  difference  t-t1 .  Thus,  assuming  the  ionosphere  to  have 
stationary  and  homogeneous  fluctuations,  we  write  the  correlation  function 


p(x-x' ,t-t')  *  <y(x,t)  m(x' ,t‘ )>  . 


(2.9) 


The  fluctuations  y(x,t)  are  further  assumed  to  have  a  power  law  spectrum,  i.e., 
for  small  Ax  =  (Ax,Ay,Az)  and  At  *  (t-t^),  we  expand  p(Ax,At)  as 


p(AX.At) 


(2.10) 


where  p,q,s  and  w  are  the  power  law  exponents  and  are  in  general  unequal.  For 

the  ionosphere  the  quantities  l_,  L  ,  L  and  t  can  vary  over  a  wide  range 

i  x  y 

depending  on  the  properties  of  the  scattering  layer  and  the  location  and  orien¬ 
tation  of  ax  with  respect  to  the  geomagnetic  equator.  A  further  restriction  on 
y(x,t)  which  is  usually  made  is  that  it  is  either  a  Gaussian  random  field  or 

that  kW>*‘  >’5  «l[with  lmax  *  max(L2,  lx,  Ly)]  in  which  case  the  distribution 

2? 

need  not  be  specified  . 

The  conditions  for  the  validity  of  the  parabolic  equation  (2.4)  can  be 

77 

expressed  in  terms  of  the  correlation  length  L(MX  and  the  correlation  time  x  . 

The  conditions  are:  (1)  that  the  spatial  scales  for  the  fluctuations  are  very 

much  greater  than  the  wavelength  of  the  source  wave  kL  ,>>  1;  (2)  kL _ «wx  ; 

max  max 

and  that  (3)  the  root  mean  square  (rms)  multiple  scattering  angle  should 
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be  small,  i.e. ,  [<y2>  R/ Lmax  ]l5«  1.  Here,  R  is  the  total  distance  of 
propagation. 

( i i )  Characterization  of  the  Scintillation 

The  scintillation  characteristics  of  a  wave  propagating  through  a  random 
medium  can  be  described  in  terms  of  single  point  and  multipoint  functions.  By 
a  single  point  function  we  mean  quantities  that  depend  on  the  distribution  of 
the  electric  field  at  a  single  point,  e.g.,  quantities  such  as  the  mean  and 
the  variance  of  the  phase  or  the  amplitude.  Multipoint  functions  such  as  corre¬ 
lation  functions  of  the  amplitude  depend  on  the  joint  probability  distribu¬ 
tion  of  the  electric  field  at  more  than  one  point.  Single  point  functions  are 
easier  to  evaluate  than  multipoint  functions  and  provide  useful  guidelines  for 
theoretical  formulation  of  approximate  solutions  of  the  problem. 

For  an  Isotropic  medium  with  a  single  scale,  Tatarski 17  has  introduced 
two  parameters  to  quantitatively  specify  the  saturation  characteristics  of 
the  propagating  wave  in  terms  of  the  medium  characteristics.  The  first  param¬ 
eter  is  the  rms  phase  fluctuation  as  computed  in  first  order  geometric  optics. 


-  k2R / p(|z|,  o)  dz  +0(Lz/R)  .  (2.11) 

•CD 

In  (2.11),  the  direction  of  propagation  has  been  taken  to  be  along  the 
z-axts.  The  parameter  $  Is  a  measure  of  the  strength  of  the  phase  fluctuations 
for  a  total  propagation  distance  R.  The  second  parameter,  called  the  diffrac¬ 
tion  parameter  (closely  related  to  Tatarski 's  wave  parameter),  is  the  square 
of  the  ratio  of  the  Fresnel  zone  size  to  the  scale  size  L  and  is  given  by 
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When  $  is  less  than  unity,  or  less  than  A"*,  either  geometric  optics  or  the 

23 

Rytov  approximation  is  valid  .  In  this  case,  there  is  a  simple  relation 
between  the  statistics  of  E(x,t)  or  of  tn  E(x,t)  and  that  of  u.  When  both  ® 
and  A*  are  greater  than  unity,  the  fluctuations  in  E(x,t)  saturate?  In  this 
case,  the  properties  of  E(x,t)  are  governed  by  statistical  considerations 
rather  than  by  the  detailed  properties  of  y.  The  path  integral  formulation  is 


particularly  useful  for  moment  calculations  in  this  regime. 

The  scintillation  characteristics  can  thus  be  summarized  in  the  AO 
diagram  as  shown  in  Figure  1.  There  are  three  regions  in  this  diagram  and 
different  theoretical  approximations  apply  in  each.  In  the  geometrical  optics 
region  the  Born  approximation  is  adequate.  When  A  is  larger  than  unity,  the 
Rytov  extension  is  adequate  as  long  as  0*<1.  In  the  saturated  region,  multiple 
scattering  effects  have  to  be  included.  A  thorough  discussion  of  the  statis¬ 
tical  properties  in  the  various  scattering  regimes  with  particular  emphasis  on 
the  relation  of  these  properties  to  the  ray  theory  can  be  found  in  Oashen's 
article22. 

(ill)  Scintillation  Characteristics  of  Ionospheric  Scattering 

Characterizing  the  scintillation  caused  by  ionospheric  irregularities 
is  complicated  because  of  the  presence  of  multiple  scales^.  Further,  the 
Irregularities  in  the  F-layer  are  highly  anisotropic®.  These  two  features  play 
a  crucial  role  In  the  overall  scintillation  behavior.  Experimental  evidence 
suggests  that  the  ionospheric  Irregularities  have  a  power  law  spectrum  with 

spectral  index  \><4.  The  inner  and  outer  scale  cutoffs  are  not  known  precisely 
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but  the  latter  Is  presumed  to  be  quite  large.  The  inner  scale  can  usually  be 

extended  to  zero  without  causing  any  major  complication.  Even  with  the  knowledge 

of  these  two  cutoffs,  the  precise  form  of  the  spectrum  (it  is  a  power  law  only 

asymptotically  for  large  spatial  frequencies)  is  not  known.  The  anisotropy 

factors  have  been  inferred  from  both  remote  and  in  situ  measurements .  In  the 

equatorial  region,  the  irregularities  are  aligned-  along  the  earth's  magnetic 

field  and  are  elongated  along  the  field  llnes^  These  are  frequently  referred 

to  as  being  rod-like  and  can  have  elongation  ratios  of  up  to  100.  In  the 

auroral  region  there  Is  considerable  evidence  that  the  irregularities  are 
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elongated  along  two  axes  and  are  sheet-like  .  These  anisotropy  factors 
play  a  major  role  In  the  observed  differences  in  the  scintillation  character¬ 
istics  of  the  polar  and  equatorial  regions. 

Two  additional  factors  that  complicate  Ionospheric  scintillation  studies 
are  the  satellite  motion  and  the  nature  of  the  data  processing.  The  collection 

and  processing  of  the  wide-band  satellite  data  is  discussed  in  a  series  of 

9-13 

reports  and  papers  by  Fremouw  et  al.  .  Two  points  regarding  the  processing 

which  are  particularly  relevant  as  far  as  this  report  is  concerned  are  that 

1)  the  raw  data  Is  subjected  to  a  sharp  cutoff  high-pass  filter  with  a  cutoff 

frequency  tc"^  and  2)  the  satellite  data  are  analyzed  in  segments  tp,  which  are 

20  to  80  seconds  long.  As  pointed  out  by  Fremouw  and  Rino,  these  two  parameters 
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play  a  decisive  role  In  determining  the  rms  phase  . 

The  phase  fluctuations  are  affected  quite  strongly  by  the  large-scale 
spatial  perturbations,  and  the  outer  scale  cutoffs  become  quite  Important. 

This  problem  does  not  occur  In  the  calculation  of  A  because  Fresnel  filtering 
effectively  suppresses  these  large  scales.  For  the  phase  fluctuations,  however, 
the  outer  scale  cutoffs  are  set  either  by  the  medium  or  by  some  measurement 
window.  In  the  wide-band  satellite  measurements,  the  scale  set  by  the 
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measurement  window,  which  is  much  smaller  than  the  actual  outer  scale  of  the 

ionosphere,  is  the  relevant  scale  cutoff.  This  scale  cutoff  can  be  taken  to  be 

some  effective  length  Lg  traversed  by  the  satellite  within  the  smaller  of  the 

two  periods  tc  or  tq.  Thus,  the  measured  phase  fluctuations  depend  on  Lg,  which 

is  a  product  of  either  tc  or  xD  and  the  effective  velocity  of  the  satellite 

with  respect  to  the  receiver.  The  effective  velocity  can  be  quite  strongly 

dependent  on  the  location  of  the  receiver  because  of  the  anisotropy  of  the 
13 

ionosphere.  Fremouw  has  argued  that  this  variation  of  the  effective  velocity 
can  explain  the  differences  in  the  scintillation  features  observed  between 
equatorial  and  polar  stations. 

The  scale  cutoff  set  by  the  measurement  window  does  not  significantly 
alter  the  amplitude  or  the  intensity  fluctuations  because  of  Fresnel  filtering 
but  the  phase  variance  is  significantly  altered.  Because  of  the  smaller  effec¬ 
tive  length  Lg,  there  is  a  significant  difference  in  the  results  for  the  phase 

variance  as  computed  in  the  Born  approximation  and  the  Rytov  approximation. 

9-13 

Rlno  and  Fremouw  compute  the  phase  In  the  Born  approximation  neglecting 
diffraction  effects.  It  is  shown  In  the  next  two  sections  that  this  can  lead 
to  a  sighificant  error  for  the  phase,  even  when  the  overall  scintillation  is 
small . 
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3.  THE  FIRST  BORN  VS.  THE  RYTOV  APPROXIMATION 

Theoretical  analysis  of  the  wlde-band  satellite  data  has  been  carried 

out  by  several  workers  assuming  that  the  scattering  Is  weak  and  that  the  Born 
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approximation  Is  adequate  In  the  computation  of  the  variance  of  the  phase 

in  the  Born  approximation,  the  propagation  effects  are  neglected.  The  result 

for<(<t>2^1s  Identical  to  that  computed  In  a  thin  phase  screen  calculation33. 

The  latter  calculation  can  be  carried  out  easily.  The  same  quantity  can  also 

be  computed  In  the  Rytov  approximation.  These  two  results  differ  considerably 

for  two  typical  satellite  passes,  over  the  equatorial  and  polar  stations, 

respectively.  To  see  the  origin  of  these  differences,  we  first  compute  the 

phase  in  the  thin  phase  screen  approximation. 

The  phase  change  induced  by  the  plasma  density  fluctuations  is  given 

33 

in  the  zeroth  order  approximation  by  the  line  Integral 

-reAyANe(r)d*  ,  (3.1) 


where  we  have  assumed  that  both  the  measurement  frequency  and  the  reference 
frequency  are  well  above  the  plasma  frequency,  rfi  Is  the  classical  electron 
radius, X  is  the  wavelength  of  the  Incident  radiation,  Afle  is  the  local  electron 
density  perturbation,  and  the  line  Integral  Is  along  a  straight  line  path  from 
the  source  to  the  receiver.  Let  the  layer  thickness  be  AL,  0  be  the  direction 
of  propagation  with  respect  to  the  z-axls  and  a^  be  the  unit  vectors  transverse 
to  the  direction  of  propagation  for  a  receiver  point  with  coordinates  (£,z) 

(See  Figs.  2,  3).  The  line  integral  (3.1)  can  be  written  as 


5$(£,z)  «  -reA  sece^  ANe(o.  +  tane  a^  (n-z),n)dv 


(3.2) 
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Then  by  direct  computation 


<£-o')  *  <^4>(e..2)<54>(e.'  .*))> 


=  r|x2  sec29AL 


AL 

x  J 0  -  )R(e.-e.'+n  tan6  ifc»n)dn 

-AL 


(3.3) 


where  R(e.,z)  Is  the  three-dimensional  autocorrelation  function  of  ANg(r). 
Introducing  the  Fourier  transform  F(Jc, kz)  of  R(fj,z)  as 


R(fi-,Z)  "  ^ifff e1("^  +  ^  Z)F(-’kz)d2|cdlcz 

In  (3.3)  and  performing  the  n  integration,  we  arrive  at  the  result: 


(3:4) 


)  *  re2x2  sec29AL 


tan9)n+  ^d-ini/AUdn  . 


(3.5) 


Assuming  that  the  correlation  length  Lz  along  the  z-di recti  on  Is  much  shorter 
then  the  layer  thickness,  we  can  neglect  the  factor  |n|/AL  and  take  the  limits 
of  Integration  to  ».  The  Integration  over  n  then  yields  the  6-functlon 

A 

27r6(kz  +  Joa^  tan  9).  Performing  the  fcz  Integration  then  leads  to  the  formulae 


R^(fir£' )  *  re2  sec29AL  ff  ^ jz  'taneik  *£)•  (3,6) 

When  the  cut offs  are  Imposed  by  the  detrender  the  Integration  over  the 
transverse  spatial  frequencies  k_  must  be  carried  out  only  over  those  spatial 
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frequencies  |kj  >  kg  ~  2tt/ Lfi .  Thus,  the  cutoffs  in  the  integration  in  (3.6) 
are  not  imposed  solely  by  the  medium  but  also  by  the  nature  of  the  data 
processing. 

We  proceed  further  by  assuming  the  same  functional  form  of  F(£,kz)  as 
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that  assumed  by  Rino  and  Fremouw  They  have  shown  that  for  a  general 
34 

anisotropy  model  (Fig.  4) 


F(k,-tane  ^  •  k)  *  ab  <AN|>Q([  Ak*  +  Bkxky  +  CkJ]‘s)  (3.7) 


where  a  and  b  are  the  axial  ratios  along  the  principal  irregularity  axes, 
respectively,  and  kx,  ky  are  x  and  y  components  of  the  transverse  wave  vector 
k,  respectively.  The  coefficients  A,  B,  C  depend  on  the  propagation  angle 
relative  to  the  principal  irregularity  axis  and  are  given  in  Eq.  41  of  Ref.  [24]. 
They  have  further  assumed  the  following  form  for  the  spectral  distribution  func¬ 
tion  Q(q) , 


Q(q) 


8ir,/j 


_2v-2 

r(v-H/2)  q0 
r(v-l)  (q£+q‘)v+V 


(3.8) 


where  qQ  Is  the  outer  scale  cutoff  wave  number  and  v  is  the  spectral  index  of 
the  power  law.  With  the  definition 

Cs  *  8w3/<^N|>q5v"1r(v+l/2)/r(v-l)  (3.9) 

we  have 


f(k,  -tan,  Sfc-k)  ■  abC,  (3.10) 

where  q2  «  Ak2  +  Bkxky  +  Ck2.  We  now  substitute  (3.10)  Into  (3.6)  with  the 

appropriate  change  of  variables  and  an  angular  Integration  and  arrive  at  the 

27h 

correlation  function 


Wy) 


*  r2  X2AL  sec9  GC, 


/qoQ(qy) 

(qi+q2)V 


(3.11) 


where 


G 


ab 

^  AC-BV4  coiF  ’ 


Jq  Is  the  Bessel  function  of  order  zero  and 

y  =  (cP2  -  Bpxpy  +  Apjj  V(AC  -  B2/^  . 


(3.12a) 


(3.12b) 


The  main  point  now  Is  the  choice  of  the  value  of  the  cutoff  that  Is  to 
be  Imposed  on  the  Integration  over  the  spectrum  of  plasma  density  fluctuations. 
We  assume  that  the  cutoff  to  be  Imposed,  as  far  as  the  outer  scale  Is  concerned. 
Is 

_  m  2ir 

qe  (3.13) 

where  the  effective  scale  length  le  is 


Ve 


in  terms  of  the  effective  velocity  V# 

C  V2  -  B V  V  +  AV21  ** 
x  x  y  y 


AC  -  B2/4 


(3.14) 


(3.15) 


and  te  =  m1n(tc,  tq).  This  assumption  is  motivated  by  the  fact  that  the 
satellite  samples  fluctuations  of  the  scale  L  *V  r  only.  Using  this 
assumption 
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1 

I 


vy>  -  w iL  sec  bscs y  t  (3-,6> 

It  should  be  pointed  out  that  this  assumption  is  almost  the  same  as  the  one 
12  13 

made  by  Rlno  and  Fremouw  .  Their  formula  involves  dropping  qQ  altogether 
in  (3.16),  a  procedure  which  is  valid  when  qQ  «  qg. 

In  ionospheric  measurements,  the  receiver  is  fixed  and  phase  fluctuations 
arising  out  of  the  satellite  motion  appear  as  temporal  fluctuations.  To  con¬ 
vert  (3.16)  to  a  formula  yielding  temporal  fluctuations  we  set  y  ■  Vgt,  and 


r5»(V>  ’  re*2  4L  “ce  sr 


M 


tr  ft 

J  (q 


0(qVfit)dq 


(q^«l2)v+H 


(3.17) 


Let  the  measured  phase  fluctuation  in  a  satellite  pass  be 

T_ 


■  4 


(t')df 


(3.18) 


An  ensemble  average  of  this  quantity  over  an  ensemble  of  statistically  iden¬ 
tical  ionospheres  is  given  by 

$2  -  ■  <y  y  .  (3.19) 

Tc  x  Tc 

The  quantity  $  can  be  obtained  from  the  correlation  function  (3.17)  by 
12  c 

setting  t«0,  i.e. , 
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00 


a  R^{0)  *  r2A2  AL  sece 


S-f-r* 

2V  <vV 


sdg _ 

jV+H 


2  2 

rgA  AL  secQ  GCS  1 


2tt 


TTO-1! 


I2v-1 )  (q§  +q2)v 


(3.20) 
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The  result  (3.20)  is  the  same  as  that  obtained  by  Rino  if  qQ  is  set  equal 

to  zero,  i.e.,  qQ  «  q0.  However,  if  qfi  is  comparable  to  qQ  the  full  expression 

for  <t  must  be  used. 

Tc 

This  result  cannot  be  compared  with  the  phase  variance  computed  in  the 

Rytov  approximation.  The  Rytov  approximation  can  be  arrived  at  in  a  variety 

of  fashions.  One  of  the  clearer  presentations  indicating  the  nature  of  the 

35  36 

approximations  involved  is  given  in  the  report  by  Call an  and  Zachariasen  *  . 

The  electric  field  is  expressed  as 

E(R)  a  EQ(R)eX(^,  (3.21) 

where  Eg(R)  Is  the  unperturbed  solution  of  the  wave  equation  at  R  with  a 
source  at  the  origin,  i.e.. 


(V2  +  q2)Eg(r)  »  <5(r),  and 

Jqlrl 

W  *  hw  ■ 

The  function  X(jr)  is  given  by 


(3.22) 


-4wr  r  j 

X(r)  -  /  d  r'E0(r-r')ANe(r')E0(r*).  (3.24) 

It  may  be  seen  from  (3.24)  that  (3.21)  can  be  obtained  simply  by  exponentiating 
the  perturbative  term  in  the  Born  approximation. 
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All  the  quantities  of  Interest  for  Ionosphere  scattering  may  be  computed 

easily  from  (3.24).  Let  us  first  split  the  function  X  into  real  and  imaginary 

*  36 
parts: 

X  *  X1  +  iX2  .  (3.25) 

Then 

E(R)  *  E0(R)eXl(-)  +  iX2(R)  .  (3.26) 


The  perturbed  phase  <|>(R)  is  then  equal  to  X2(R).  Thus, 

<♦*>■<*&>> 

*  -  r'<x2>]  '  (3-27) 

Similarly,  the  log-intensity  L  can  be  defined 


I  =  log  |E(R)|2-<t>  +  2X-| (R) , 
where  log  |EQ(R)|2,  and 

•  2|^X|2^Re^X2\>j. 


(3.28) 


(3.29) 


Thus,  to  evaluate  either  the  phase  variance  or  the  log-intensity  variance  we 
need  to  compute  only  ^|X|2^and  <^X2^.  These  calculations  are  carried  out 

in  Appendix  A.  The  result  is 

RM 

<^|X(R)|^>»  rj  Xzy  dx  J  dlki  F(0,kx)  (3.30) 

R' 

R" 

<X2(R)>-  -  rj  X2  J  ixj  d2kx  F(0,kx)exp(iki  /qA)  ,  (3.31) 

R' 
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where  A  =  R/x(R-x)  is  the  wave  curvature,  R'-RM  is  the  distance  traversed  in 
the  layer,  and  kx  is  the  component  of  the  spatial  wave  number  perpendicular 
to  the  direction  of  transmission. 

A  detailed  discussion  of  the  nature  of  ^|X(R)|^>  and  ^X2(R)^>  is 
given  in  Chapter  10  of  Ref.  [23].  The  quantity  <^|X(R)|2^>in  (3. 3C)  is  The 
same  result  as  obtained  by  computing  X(R)  in  geometrical  optics  (cf.  3.26). 

In  fact,  ^  I X (JR )  I2  y  is  the  same  as  &  in  geometrical  optics.  The  quantity 
<X2(R)>  in  (3.31)  is,  however,  different  from  the  value  obtained  from  a 
geometrical -optics  calculation.  The  geometrical -optics  result  for  ^X2(R)]> 
is  obtained  by  expanding  the  exponential  in  (3.31)  and  replacing  it  by  unity. 
This  would  correspond  to  the  condition  k  /qA«l  for  all  spatial  wave  numbers 
in  the  integration  in  (3.31).  Let  some  typical  kj.  be  equal  to  L  where  L 
Is  the  correlation  length.  Thus  the  condition  for  the  validity  of  geometrical 
optics  becomes  l/L2qA«l.  Since  A"1  -v  R  this  condition  is  the  same  as  the 
Fresnel  condition  A  =  (  /R/q/L)2«l .  In  terms  of  (3.30)  and  (3.31),  the  phase 
variance  is  given  by 
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The  geometric-optics  result  for^2/*  corresponds  to  replacing  the  cosine  term 

in  brackets  by  unity.  This  approximation  is  valid  if  and  only  if  the  argument 

of  the  cosine  can  be  neglected  for  all_  relevant  values  of  k  and  k  .  To  exam- 

ine  the  result  (3.33)  further,  we  replace  the  wave  curvature  function  by  its 

value  A  in  the  middle  of  the  scattering  layer.  (For  a  scattering  layer  at 

300  km  from  the  receiver  and  a  layer  thickness  of  50  km  and  a  satellite  height 

•3  —1  -3  —1 

of  1000  km,  A  varies  between  4.5  x  10  km  and  5.0  x  10  km  and  may  be 
approximated  by  some  average  value  without  introducing  major  errors.)  In  this 
approximation, 

<^<p2)>a  r2X2sec29£L J"  dkxdkyF(kx  ky,  -tan9cos<J>kx  -tan8sin<|>ky) 

i  1  k2+k2nan28(kvcos<H’M1n(<,n 

x  l  [1+cos  \-2~Z - * - * - I  ].  (3.34) 

ql 


We  are  now  in  a  position  to  incorporate  the  anisotropy  of  the  ionosphere  into 
our  expressions, 

33  34 

Using  the  generalized  anisotropy  model  of  Rino  and  Fremouw  ’  we  may 
rewrite  (3.34)  as 


2  2  2 

r“X  sec  ab 
e 


“l  ■  TW,2  f  ♦  Vy  * 

k^+k^+tan2  9(ky  cos$  +  ky  sin#)2 


x  ^£1+  cos{ 


x  y 


qA 


•>J 


(3.35) 


(3.36) 


where  A,  B  and  C  are  functions  of  the  anisotropy  scale  factors  a,b  and  the 
anisotropy  angles  (see  Fig.  4).  Eq.  (3.36)  can  be  simplified  by  a  diagonaliz¬ 
ing  transformation  followed  by  integration  over  the  angular  variables.  The 
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details  are  given  in  Appendix  B.  The  result  is 


«  — — 1 - f dkkQ(k)[l  +  cos  (k2Y] )0Q(k2Y2 )], 

4^AC  -  B2/^ 


(3.37) 


where  Yl  and  y2  are  denned  in  Appendix  B.  The  limits  of  integration  over 
the  spatial  frequencies  must  be  determined  once  again  by  the  nature  of  the 
data  processing.  As  per  our  earlier  discussion,  we  set  the  outer  scale  to 
qe  and  the  inner  scale  to  «  to  obtain 


00 

— 1—  /  dkkQ(k)[l  +  cos(k2r1)J0(k2r2)]. 
yAC  -  B2/4 j 


(3.38) 


Thus  the  phase  variance  ^|>2^> 


in  the  Rytov  approximation  is  given  by 


reX2sec2  8AL  ab  <AN 2)>  J  dkkQ(k)[l  +  cos(k2Y. ) JQ(k2Y?)]. 
4ir(AC  -  B2/4)Ss  <e  1  0  2  (3.39) 


Substituting  the  expression  Q,k)  from  (3.8)  and  scaling  out  qfi,  we  obtain 
00 

<*2>*  <3-40> 


where 


8  = 


r2X2  At  sec20  ab  <JlN 2^>  q 
4tt  (AC  -  lZH)h 


(3.41) 
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Similarly,  using  ;3.29),  (3.30),  and  (3.31),  the  variance  of  log-intensity  is 
given  by 

”  ^y2  =  28  ^(q2yqk2'k+  k2jv+»s|1"C0S^k2qeYl  ^O^eVj  (3,42) 

Note  the  relative  sign  in  the  bracketed  terms  in  (3.40)  and  (3.42). 
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4.  DATA  ANALYSIS 

Ionospheric  scintillation  data  obtained  from  satellite  observations 

have  been  analyzed  by  several  workers.  Of  particular  Interest  to  this 

6  7  12  13 

section  Is  the  work  of  Crane  *  ,  Rino  and  Fremouw  .  Crane  has  presented 
an  extensive  discussion  of  the  data  gathered  from  the  U.S.  Navy  Navigational 
System  Satellites  (NNSS)  at  Millstone.  Details  of  the  receiving  system  at 
the  Millstone  Radar  Facility  Is  described  in  his  review  article^.  Two  phase 
coherent  sources  at  150  and  400  MHz  were  used.  This  enabled  the  elimination 
of  phase  shifts  introduced  due  to  the  satellite  motion.  Crane  analyses  of 
the  phase  and  amplitude  scintillation  through  the  evaluation  of  such  measures 
as  the  phase  variance  log-intensity  variance  (*4^(^  = 

(^Tir)  a  x)  and  the  intensity  scintillation  S^.  Rino  has  carried  out  a  similar 
analysis  for  the  wide  band  P76-5  satellite  data.  The  accumulation  and  pro¬ 
cessing  of  the  wide  band  satellite  data  Is  described  In  the  paper  of  Fremouw 

Q 

et  al.  The  wide  band  satellite  used  coherent  signals  in  the  VHF,  UHF,  the 
L-Band  and  the  S-Band.  Observations  were  made  In  Ancon  (Peru),  Kwajalein 
(Marshall  Islands),  Poker  Flat  (Alaska)  and  Stanford  (California).  Rino  and 
Fremouw  present  detailed  analysis  of  and  for  each  of  these  observations 
stations. 

(1 )  Weak  vs  Strong  Scatter  Theories. 

The  main  point  we  wish  to  discuss  Is  the  theoretical  models  used  by  these 
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researchers  to  analyze  the  data.  Rino  ’  concludes  in  his  1979  papers 
that  a  major  part,  if  not  the  whole,  of  Ionospheric  scattering  phenomena 
might  be  described  by  the  Born  approximation.  His  conclusion  rests  on  the 
following  two  arguments.  First,  in  the  case  of  the  ionosphere  the  outer 
scale  is  sufficiently  large  that  ,0[xM  ]\  -here  a'1  •  15  the  reduced 
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height,  is  much  smaller  than  unity.  Second,  for  the  ionospheric  scattering, 

Rumsey's  intensity  randomization  parameter  U, which  is  essentially  the 

average  square  intensity,  is  much  smaller  than  unity.  It  is  argued  that  these 

two  factors  are  sufficient  to  justify  the  application  of  the  Born  approxima- 

13 

tion  theory  to  data.  Fremouw  has  also  given  support  to  this  argument  in 
his  1980  paper  where  he  concluded  that  U-0.1  for  ionospheric  scattering. 

Crane's  analysis^  of  the  Millstone  data  points  to  a  somewhat  different 
conclusion.  Crane  compares  the  amplitude  scintillation  data  to  both  Born 
and  Rytov  approximation  results.  In  the  Rytov  calculation,  the  appropriate 
measure  Is  <X1  >,  (cf.  3.29)  proportional  to  the  variance  of  the  log-inten¬ 
sity.  In  the  Born  approximation,  he  uses  the  intensity  scintillation  S4.  His 
comparison  of  theory  and  data  indicates  that  neither  Born  nor  Rytov  agrees 
with  the  data  in  its  entirety.  Simulatenous  observations  of  and  S4 
clearly  shows  saturation  behavior  beyond  a ^  -|50  =  5.56  dB  and  S4  15Q  values 
above  0-8.  Thus  there  does  appear  In  ionospheric  scattering  data, points 
that  cannot  be  fitted  by  either  Rytov  or  Bom  approximation.  Even  Fremouw's 
S4-  scatter  diagram  (see  fig.  11)  indicates  such  points  with  large 

values  of  S4  and  clearly  a  weak  scatter  theory  is  not  likely  to  be  correct 

12 

for  these  data  points.  The  fit  of  Rino  for  passes  such  as  Ancon  32-06  and 
Poker  Flat  12-24  show  deviations  when  S4  values  are  large  particularly 
around  the  edges  of  the  pass.  Also,  the  scintillation  observed  at  Poker 
Flat  Invariably  shows  enhanced  scintillation  at  points  In  and  around  the 
middle  of  the  pass.  This  enhancement  is  distinct  from  the  gradual  enhance¬ 
ment  observed  at  the  edges,  see  for  example,  for  the  Poker  Flat  Pass  12-24 
(Fig.  5),  Pass  (6-14)  (Fig.  6)  and  our  theoretical  Rytov  curves  for  Pass 
6-39  (Figs.  17,  18).  There  are  portions  of  the  pass  for  which  S4  orx^ 
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clearly  exceeds  the  limits  of  the  weak  scatter  theory  and  saturation  sets  In. 
Hence,  the  weak  scatter  predictions  are  much  larger  than  the  data.  In  such 
regions  there  is  no  alternative  except  to  use  strong  scatter  theory. 

(ii )  Born  vs.  Rytov  Approximation 

Even  in  regions  where  is  small  the  question  arises  as  to  whether  the 
Born  orthe  Rytov  approximation  is  to  be  used.  Crane7  compares  both  approxima¬ 
tions  to  data  and  favours  the  Rytov  approximation  (see  Fig.  7).  The  Born 
approximation  results  for  ranging  from  0.2  to  0.4  lie  below  the  data 

and  those  ranging  from  0*8  to  1  lie  above  the  data.  The  scatter  of  the  data 
around  theory  is  above  the  one  dB  limit  whereas  the  scatter  is  less  in  the 
Rytov  fit.  It  appears  from  this  analysis  that  the  Rytov  fit  is  better. 

The  phase  variance  calculations  that  we  have  carried  out  also  raises  some 
questions  about  the  Born  approximation.  The  Bom  approximation  result  for 
the  phase  as  derived  by  Rino  is  the  same  as  the  geometric  optics  results 
given  by 


<<02>b  -  <!x|2> 


(4.1) 


where  <(|X|^  is  given  by  (3.3). 

The  Rytov  approximation  result  is  given  by 

<*2>  -  i/z[<txi2>  -  rs<x2>]  (4  2) 

where  <X  >  is  given  by  (3.31).  These  two  results  are  equal  if  the  outer 
scale  is  so  large  that  for  a  significant  part  of  the  integration  in  (3.31) 
k^  is  essentially  zero,  i.e.,  diffraction  effects  due  to  the  large  scales 
are  negligible  and  dominate  the  contributions  to<$2>  .  In  ionospheric 
scattering  the  outer  scale  is  not  the  relevant  scale.  As  explained  in 
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Sec.  2,  limits  of  integration  of  the  spectral  function  in  (3.30)  and  (3.31) 


are  set  by  the  detrender  and  is  some  scale  qfi  which  may  be  much  smaller  than 
the  outer  scale  cut-off  qg.  Thus  in  the  integration  in  (3.31)  it  may  not  be 
correct  to  approximate  the  exponential  by  unity.  To  examine  these  differ¬ 
ences  closely,  we  have  computed  the  ratio  of  the  phase  variance  in  the  Rytov 
approximation  to  that  in  the  Born  approximation.  From  (3.39)  this  ratio  can 


be  seen  to  be 


<A 


^dk  k  Q(k)  1  +  cos  (k2^)  JQ(k2Y2) 

00' 

2  f  dk  k  Q(k) 


(4.3) 


The  computation  has  been  carried  out  for  the  Kwajalein  Pass  22-13  and  the 
Poker  Flat  Pass  6-39  (See  Figs.  8,  9,  10).  A  rather  interesting  conclusion 
emerges  from  this  comparison.  For  the  Poker  Flat  case  there  is  very  little 
difference  between  the  Rytov  and  Born  results  for  sheet-like  irregularity 
ratios  of  10:10:1  and  40:40:1.  However,  there  is  a  maximum  of  10*  differ¬ 
ence  for  irregularities  with  ratios  40:10:1  l.e. ,  those  that  are  Intermediate 
between  rods  and  sheets.  Thus  for  the  usually  accepted  sheet-like  model  for 
the  irregularities  the  Bom  approximation  is  close  to  the  Rytov  results. 


In  the  Kwajalein  case  irregularities  are  rod-like  and  there  is  significant 
difference  between  Bom  and  Rytov  results.  For  an  Irregularity  ratio 
10:1:1  there  Is  a  maximum  deviation  of  about  15*  which  appears  quite  sharply 
within  the  satellite  pass  (see  Fig.  10).  For  higher  Irregularity  ratios 
of  40:1:1  and  100:1:1  the  deviations  become  more  pronounced.  The  100:1:1 
result  lies  slightly  below  40:1:1  but  has  a  very  similar  shape.  All  these 
results  point  to  the  conclusion  that  there  is  significant  difference  between 
Bom  and  Rytov  results  in  the  equatorial  stations  which  again  imply  that  the 


diffraction  effects  may  not  be  negligible.  A  reason  why  these  effects  are 
more  pronounced  In  the  equatorial  station  compared  to  the  polar  station  Is 
the  rod-like  nature  of  the  irregularities  in  the  former  compared  to  the 
sheet-like  Irregularity  In  the  latter.  For  the  two  passes  we  considered 
(and  also  for  other  typical  passes  of  the  P  76-5  satellite)  the  x-component 
of  the  satellite  velocity  (in  the  geomagnetic  north-south  direction)  is 
ten  times  larger  than  the  y-  component  (in  the  east-west  direction).  When 
the  irregularities  are  sheet-like  the  effective  velocity  Vg  is  still  dominat 
ed  by  the  Vx  component  and  is  quite  large.  The  effective  outer  scale 
qfi  ■  Lg1  *(veTe)-1  (See  Sec.  2)  is  then  quite  small  and  diffraction  effects 
are  small.  When  the  irregularities  are  rod-like,  the  effective  velocity  is 
considerably  reduced  (it  is  only  about  15%  of  Vx).  In  this  case  the  effect¬ 
ive  outer  scale  Vfi  is  large  and  the  diffraction  effects  are  non-negliglble. 

There  is  experimental  evidence  to  indicate  that  diffraction  effects 
are  lowering  the  phase  variance  in  equatorial  stations.  Fremouw  examines 
the  S^-o^  diagram  for  several  passes  (see  Fig.  11).  He  finds  that  the  ratio 
S^-cr^  is  consistently  higher  in  equatorial  stations  as  opposed  to  polar 
stations.  Fremouw  seeks  an  explanation  for  this  effect  in  such  factors  as 
the  incidence  angle  0,  ionospheric  anisotrophy,  source  motion  etc.  He 
concludes  that  it  is  the  course  motion  coupled  with  the  effect  of  detrending 
on  phase  scintillation  that  explains  this  difference  in  the  ratio 
Our  calcuatlons  would  not  only  support  this  view  but  also  strengthen  it. 
Fremouw1 s  viewpoint  is  that  there  is  more  diffraction  introduced  because  of 
the  source  motion  and  the  detrender  cutoffs  in  the  equatorial  case  compared 
to  the  polar  case.  He  has,  however,  examined  the  diffraction  effects  in  S4 
only  and  not  included  it  in  his  calculation  which  is  essentially  a 
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geometric-optics  calculation.  When  we  carry  out  the  calculation  in  the  Rytov 

approximation  o  decreases  because  of  diffraction.  Thus  the  increase  of  the 
♦ 

ratio  $4/0^  arises  not  only  because  of  the  increase  of  S4  but  also  because 
of  the  decrease  of  due  to  diffraction.  The  decrease  in  as  pointed  out 
earl ier, appears  only  in  the  Rytov  approximation. 

( i i i )  Rytov  Calculations 

To  examine  the  sensitivity  of  the  phase  and  intensity  scintillation 
we  have  computed  these  quantities  in  the  Rytov  approximation  for  the  Kwajalein 

Pass  22-13  and  the  Poker  Flat  Pass  6-39.  For  this  purpose  we  have  arbitrar- 

2 

ily  set  the  value  of  <$  >  to  0.2  in  the  middle  of  the  pass.  The  variance  of 
log-intensity  is  computed  with  the  phase  variance  normalized  to  this  value. 

We  have  chosen  two  values  of  the  spectral  index  v*  1.3  andv*  1.5  which  are 
representative  of  the  values  encountered  in  the  Ionospheric  measurements. 

We  first  discuss  the  equatorial  case.  We  have  chosen  three  anlsotrophy 
ratios  10:1:1.  40:1:1  and  100:1:1.  The  phase  variance  <$>  monotoni call y 
Increases  for  Increasing  anlsotrophy.  (See  Figs.  12  &  13).  It  also  increases 
smoothly  towards  the  edges  of  the  pass,  an  effect  caused  primarily  by  the 
sec  9  factor  in  (3.41).  The  calculated  phase  variance  Is  not  likely  to  be 
correct  towards  the  edge  of  the  pass. if  the  corresponding  log-intensity 
variance  is  high.  Then  the  phase  and  the  intensity  should  be  given  by  the 
strong  scatter  formulae.  Crane  finds  that  experimentally  this  occurs  at 
<1  >  -  <\>  - 1  - 2.  Since  our  normalization  is  arbitrary  we  cannot  find 

from  our  graphs  at  what  point  the  weak  scatter  approximations  breakdown. 

Our  calculations  show  that  it  is  possible  to  distinguish  anistropy  ratio 
of  10:1:1  from  40:1:1  but  the  latter  anisotropy  yields  results  quite  close 
to  that  of  100:1:1.  Thus  It  is  difficult  to  find  accurate  estimates  of  the 
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anisotropy  ratios  from  <♦  .>  and  <i  >-<t>  .  These  results  are  also  quite 
Insensitive  to  the  spectral  indices. 

The  Poker  Flat  pass  exhibits  a  more  complicated  behavior.  Vie  have 
analyzed  three  anisotropy  ratios  10:10:1,  40:10:1  and  40:40:1.  The  variance 
<*2>  is  once  again  arbitrarily  set  to  0.2  slightly  to  the  left  of  the  middle 
of  the  pass.  The  most  striking  feature  is  the  abrupt  increase  of  the  phase 
variance  in  the  middle  of  the  pass  for  all  the  anisotrophy  ratios.  This 
variation  is  in  fact  spurious  and  is  caused  by  a  very  rapid  variation  in 
the  azimuth.  At  these  points  the  error  bars  on  the  measured  velocities 
are  very  large.  So,  the  computation  of  these  data  points  is  far  from 
accurate.  These  data  points  appear  at  the  top  of  each  curve  (see  Fig.  16-18). 
Excluding  these  data  points  still  leaves  a  jagged  curve  representing  the 
phase  and  the  intensity  data.  The  Intensity  scintillation  shows  quite  a  bit 
of  variation  and  the  weak  scatter  theory  probably  does  not  apply  at  the 
scintillation  peaks.  In  any  case,  the  sharp  Increases  observed  in  the 
intensity  data  (see  Fig.  6)  and  our  theoretical  curves  (Fig.  17,  18)  argue 
against  the  exclusive  use  of  the  weak  scatter  theory. 
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CONCLUSION 

We  sunmarize  the  major  conclusions  as  follows. 

12  13 

1.  Contrary  to  the  claim  of  Rino  and  Fremouw  the  application 
of  the  weak  scattering  theory  to  the  entire  ionospheric  data 
base  is  not  justified.  In  both  equatorial  and  polar  stations, 
there  appear  data  points  for  which  the  variance  of  the  inten¬ 
sity  is  quite  large.  In  polar  stations  such  points  appear 
when  the  signal  trajectory  lies  transverse  to  the  short  axis. 

In  equatorial  stations  high  values  of  are  encountered  at 
the  edges  of  the  pass.  Under  such  conditions  a  strong  scatter 
theory  should  be  used. 

2.  Even  in  regions  where  weak  scatter  theory  is  applicable,  the 
Rytov  approximation  is  preferable  to  the  .Born  approximation 
for  the  phase  calculations.  This  observation  is  consistent 
with  Crane's  analysis7  of  the  Millstone  data.  In  equatorial 
regions  there  is  a  very  significant  difference  between  the 
Rytov  and  Bom  approximation  results  indicating  that  the 
diffraction  effects  are  non-negligible. 

3.  The  diffraction  effects  reduce  the  phase  variance  and  lead  to 
smaller  cr.,  an  effect  which  is  particularly  pronounced  for 
equatorial  stations.  The  Itytov  approximation  partially  accounts 
for  these  diffraction  effects  whereas  the  usual  geometric  optics 
calculation  neglects  them  altogether.  Thus  for  the  same  strength 
of  electron  density  perturbations,  S^/o^  is  larger  in  equa¬ 
torial  than  in  polar  stations.  This  is  another  theoretical 
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reason  that  can  be  advanced  to  explain  the  enhancement  in 
observed  by  Fremouw13  in  equatorial  stations  compared  to  that  in 
polar  stations. 
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Appendix  A 
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(A. 3) 


^here  (x^ ,  x.,  xk)  are  the  components  of  x  and  yk  has  been  neglected.  Thus 
the  exponent  in  (A.l)  reduces  to  exp  iq{(x.y.  +  x -y  ■)-— ,-D— ---r}. 

1  *  j  j  aj^vk  •  Xj^y 


In  terms  of  these  new  variables  (A.l)  can  be  rewritten  as 

,2 


<|X(R)|2>  =  r2 


U  /d2x - - 

J  J  (R_x» 


(R-xk)2x2 


xexp  h(xiyi+x.y.)7-T|rj-r| 


(A. 4) 


An  integration  over  dx.  dx.  yields 
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(A.5) 


Substituting  the  result  (A.5)  in  (A.4),  we  arrive  at 


R' 

<|X(R)|2>  =  r2  X2  /  dx  f  dy  t  RglyJ.  (A. 6) 

R" 

where  y  is  the  integration  variable  parallel  to  the  straight-line  path  and 

II 

(R1  -  R")  is  the  distance  traversed  in  the  layer.  In  terms  of  the  frequency 

spectrum  of  can  rewrite  (A. 6)  as 

R' 

<|X(R)|2>  *  r2  \2  J  dx  /d\  F(0,kA)  (A. 7) 

0 
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An  entirely  analogous  stationary-phase  calculation  can  be  carried  out 
o 

for  evaluating  <X  (R)>.  From  (3.7)  and  (3.8)  we  obtain 

f  |  R|  2 

<*2(*)>  •  rg  /a3 V"  d3-"|R.  r'llf-  r“||r‘||r“|R(^  ' 

x  exp  iq{  |R  -  r'  |  +  |jr*  |  +  | R  -  r“  |  +  |  r"  |  -  2| R| } .  (A. 8) 

Once  again,  the  stationary  path  is  the  straight  line  from  the  origin  to  the 
point  Jl.  Expanding  the  exponential  as  before, 

| R  -  r' |  +  | r' |  +  | R  -  r"|  +  |r“|  -  2|R| 


■  ««? ♦  »?/*)  *  (*J *  ^4»  * 


(A. 8)  can  be  rewritten 


9  ?  /  dx.  R 

<X^(R)>  =  1 


R"  J  J 

exp[iq{x^  *  y3/4  *  x3  *  y3/4)(i  * 


(A. 9) 


Performing  the  Gaussian  integration  ,  we  obtain 

2  r  r 

9  ir>  /  dx.R  /  - 

<X2(R)>  =  — /  (ftk,~k)Xk  J  d2  yA  dy(iR(yi(,^) 


y2  R" 
expCic^K^  +  _I_}3, 


(A. 10) 


(A. 11) 


2  2  2 

where  yx  *  y^  +  y^.  Introducing  the  Fourier  transform  of  R(r)  as  before, 
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<X  (R)>  - - 2 

q 


R' 

4-2r-&-  ^x/d^FCO.kJexpCiqkfc^^-)], 
R" 


where  the  subscript  k  on  x^  has  been  dropped. 


! 
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Appendix  B 


Consider  the  following  transformation. 


(B.l) 


which  diagonalizes  the  quadratic  form  of  the  argument  of  Q  in  (3.17).  In 

terms  of  the  variables  k',  k'  the  integral  over  the  second  term  in  (3.15)  =  I 

x  y 

is 
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k"  -ySkT 


=  k  cos  e 


c»  =\[i- 


BV4A  k^ 


k"  s  k  sin  0 


successfully,  we  obtain 


where 


!  =  _L  (AC  -  B2/4)_is  tJk"k**Q(k") 

2ir 

x ^decos  ^(y-j+Y2C0S  ^  20~<|>) )  k j  . 

0 

% 

T’*iP+rfe) 

^sip-rk) 

«2  AV^ 


tan'1# 


C  -  BV4A 


ai(C  -  B  /4A)  -  a3A 
Performing  the  integration  over  0,  we  arrive  at 


I  * 


\J  AC-B^/4 


/ 


cos(k2Yl }J0(kZY2)Q(k)kdk, 


(B.4) 


(B.5) 


(B.  7) 
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FIGURE  4:  Ionospheric  Anisotropy  Geometry:  s,  r,  t  define  an  orthogonal 
coordinate  system,  s  Is  along  the  magnetic  field  at  an  angle 
i p  wr»t  the  x-axIs.  The  cross  hatched  region  defines  the  plane 
of  magnetic  meridian.  The  other  elongation  axis  r  Is  at  an 
angle  6  wr«t  the  plane  of  the  magnetic  meridian. 
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Data  from  Poker  Flat  pass  12-24  showing  Isolated  geo¬ 
metrical  enhancement  attributed  to  sheetlike  irregu¬ 
larities.  (reproduced  from  C.  L.  Rino,  Rad.  Scl.  1_4* 
1135  (1979)). 
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t«h.  in 
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FIGURE  6:  Data  from  Poker  Flat  pass  6-14  together  with  theoretical 
curves  showing  evidence  of  height  change,  (reproduced 
from  C.  1.  Rino,  Rad.  Sci.  14j_  1135  (1979)). 


Univ.  Time 

FIGURE  8:  Ratio  of  phase  variance  computed  In  the  Rytov  and  Born  approximations  for  Poker  Flat  Pass  6-39 
with  v*l .3  and  asymmetry  ratios  ***.10:10:1,  ...40:10:1  and— 40:40:1. 


FIGURE  1 


:  Scatter  diagram  of  scintillation  indices  for  intensity 
(S.)  and  phase  (ad>)  for  50  VHF  data  segments  ranging 
inHength  from  20  to  85s.  Note  consistency  of  pattern 
from  two  equatorial  stations,  Kwajalein  and  Ancon,  and 
separate  pattern  formed  by  data  from  a  high-latitude 
station.  Poker  Flat,  (reproduced  from  E.  J.  Fremouw, 

J.  Atmos.  Res.  42^  775  (1980)). 
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UM1U.  T1IC 

Mean  square  phase  vs.  univ.  time  for  Kwajalein  t>ass  22-13  with  v=1.3  and 
asymmetry  ratios  ***10:1:1,  ...40:1:1  and  — 100:1:1. 


UMIU.  Tint 

FIGURE  14:  Variance  of  log  Intensity  vs.  unlv.  time  for  Kwajaleln  Pass  22-13  with  v=1.3 
and  asymmetry  ratios  ***10:1:1,  ...40:1:1  and  — 100:1:1. 


Mean  square  phase  vs.  univ.  time  for  Poker  Flat  Pass  6-39  with  \*=1.5 
and  asymmetry  ratios  ***10:10:1,  ...40:10:1  and  — 40:40:1. 


UNXU.  TINE 

FIGURE  18:  Variance  of  log  of  Intensity  vs.  univ.  time  for  Poker  Flat  Pass  6-39  with 
v=1.5  and  asymmetry  ratios  ***10:10:1,  ...40:10:1  and  — 40:40:1 


